1+1 


National  Defense 
Defence  natlonale 


DTIC 


Fl  _'CTE 
JUL  2  4 1991 


AD-A239  039 

!  Ill  Hi  ilia  i,...  .....  . 


AN  ANALYSIS  OF  COMPLEX  DOWN  MIXING, 
FILTERING  AND  DECIMATION  IN  ACES  (U) 


by 

Eric  April 


SKTEisirnoi*  sTAnagsT  k  ~\ 

Appiovad  for  pus uc  rsieesat  1 
Clalirc-ttsd  ; 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 

TECHNICAL  NOTE  91-2 


Canad'ai 

y  i  v 


91-05993 


i 

I 


December  1990 
Ottawa 


1*1 


National  Defense 

Defence  naltonale 


AN  ANALYSIS  OF  COMPLEX  DOWN  MIXING, 
FILTERING  AND  DECIMATION  IN  ACES  (U) 

by 

Eric  April 

Communications  Electronic  Warfare  Section 
Electronic  Warfare  Division 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 

TECHNICAL  NOTE  91-2 


PCN 

041LK11 


December  1990 
Ottawa 


ABSTRACT/RESUME 


ABSTRACT 

This  report  deals  with  the  problems  related  to  the  use  of  an  all-digital  HF/VHF 
wideband  receiver  in  a  digital  communication  intercept  system.  The  high  output  rates 
and  large  bandwidths  associated  with  these  receivers  result  in  the  need  for  a  post¬ 
processing  system  which  can  reduce  the  sampling  rate  and  isolate,  in  frequency,  the 
bandwidth  of  interest.  Therefore,  a  means  of  digitally  filtering  and  mixing  the  signal  to 
baseband  is  required.  Additionally,  a  decimation  process  should  be  included  in  the 
system.  General  and  detailed  descriptions  of  the  Mixing,  Filtering  and  Decimation 
System  (MFDS)  are  given.  Simulations  of  the  MFDS  are  also  provided. 


RESUME 

Ce  rapport  traite  des  problemes  relies  a  1'utilisation  de  recepteurs  a  large  bande 
comptetement  numeriques  fonctionnant  dans  les  bandes  HF  et  THF  a  l'interieur  d'un 
systeme  d'interception  de  communications  numeriques.  Les  taux  d'echantillonnage  de 
sortie  excessifs  ainsi  que  les  largeurs  de  bande  trop  etendues  associes  a  de  tels  recepteurs 
conduisent  de  faqon  imminente  au  besoin  de  conceptualiser  un  systeme  rattache  a  la 
sortie  du  recepteur  qui  reduira  le  taux  d'echantillonnage  et  isolera,  en  frequence,  la 
bande  spectrale  d'interet.  Done,  des  methodes  numeriques  de  filtrage  et  de  mixage  du 
signal  &  la  bande  de  base  sont  requises.  De  plus,  une  operation  de  decimation  devra  etre 
incorporee  au  systeme.  Une  description  generate  puis  detaillee  du  Systeme  de  Mixage, 
de  Filtrage  et  de  Decimation  (SMFD)  sont  donnees.  La  validite  du  systeme  (SMFD)  est 
evaluee  suite  aux  simulations  par  ordinateur. 
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EXECUTIVE  SUMMARY 


This  report  elaborates  a  system  which  reduces  the  sampling  rate  of  a  signal 
coming  from  a  digital  wideband  receiver  (which  has  ~  4  MHz  bandwidth)  and  isolates, 
in  frequency,  the  bandwidth  of  interest.  This  system  is  called  the  Mixing,  Filtering  and 
Decimation  System  (MFDS).  This  work  was  carried  out  in  support  of  research  into  the 
Advanced  Communications  Electronic  Support  Measures  System  (ACES). 

The  choice  for  the  kind  of  digital  filters  to  be  used  in  the  MFDS  is  dependent  on 
three  specifications:  (1)  the  transition  width  defined  as  the  spectral  region  going  from  the 
passband  to  the  stopband,  (2)  the  stopband  attenuation  and  (3)  preservation  of  the  phase 
linearity  of  the  filter.  From  these  specifications,  an  FIR  (Finite  impulse  Response)  filter 
using  a  Kaiser  window  weighting  function  with  255  taps  was  selected.  Detailed 
specifications  for  this  filter  are  contained  in  the  report. 

Some  typical  values  associated  with  the  receiver  has  to  be  chosen.  These  include 
the  receiving  bandwidth  (4  MHz),  sampling  rate  (8  MHz)  and  scanning  range  capability 
(=  3-300  MHz  ie.,  HF/VHF  bands).  From  there,  an  analysis  of  the  overall  system  has  led 
to  the  need  (mainly  due  to  filter  limitations)  for  two  decimation  operations  such  that  a 
new  lower  sampling  frequency  is  defined.  Following  this,  a  detailed  description  of  the 
MFDS  showed  three  things:  (1)  how  to  build  a  complex  bandpass  filter,  (2)  how  to 
perform  a  complex  down-mixing  using  complex  bandpass  filters  and  finally,  (3)  under 
which  conditions  a  decimation  operation  has  to  be  performed  in  order  to  prevent 
aliasing. 

Hardware  implementation  issues  are  also  considered.  One  important  issue  is  the 
need  for  very  powerful  processors  to  meet  the  real-time  computational  demand.  Since 
common  DSP  chips  (e.g.,  TMS320C30)  do  not  manifest  such  power,  INTEL  iYVarp 
parallel  processor  boards  are  proposed  as  one  of  the  best  solutions.  In  addition,  because 
the  signal  of  interest  is  not  completely  isolated  from  other  signals  even  after  two  levels 
of  filtering  as  noted  above,  further  processing  will  be  required.  This  latter  issue  will  be 
discussed  in  a  subsequent  report. 

Finally,  computer  simulations  have  been  used  to  prove  the  validity  of  the  theorv 
used  under  the  MFDS. 
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1.0  INTRODUCTION 

ACES  is  the  acronym  for  Advanced  Communications  Electronic  Support  Measures 
System.  It  is  a  system  intended  to  integrate  as  much  as  Electronic  Support  Measures 
(ESM)  functions.  These  include  signal  interception,  signal  detection,  direction  finding, 
modulation  recognition,  demodulation,  etc...  This  system  will  utilize  all-digital  state-of- 
the-art  technologies  such  as  wideband  digital  receivers  and  high-speed  digital  signal 
processors.  For  the  sake  of  the  present  work,  a  sampling  frequency  of  8  MHz  and  a 
receiving  bandwidth  of  4  MHz  will  be  assumed  for  the  receiver.  It  will  also  be 
presumed  that  this  digital  receiver  will  be  able  to  sweep  the  received  bandwidth  through 
the  whole  HF/VHF  band. 

Trying  to  process  signals  at  a  sampling  frequency  fs  of  8  MHz  is  a  really  difficult 
task  in  a  real-time  digital  communication  intercept  system.  In  fact,  looking  at  a  4  MHz 
bandwidth  (BW=fs/2)  inside  the  HF/VHF  bands  leads  to  the  observation  of  more  than 
one  signal  of  interest  (SOI).  Assuming  the  following  types  of  modulation, 

a) HF  band  (3  -  30  MHz):  (with  3  KHz  channel  spacing) 

Upper  Side  Band  (USB)  :  2.7  KHz  bandwidth 
Lower  Side  Band  (LSB)  :  2.7  KHz  bandwidth 
On-Off  Keying  (OOK)/  Amplitude  Shift  Keying  (ASK)  :  30  to 
50  Bauds  -  150  Hz  bandwidth 

Frequency-Shift  Keying  (FSK)  :  45/100  Bauds,  250/300  Hz 
bandwidth 

Quadrature  Phase-Shift  Keying  (QPSK),  Binary  Phase-Shift 
Keying  (BPSK) 

Double  Side  Band  -  Suppressed  Carrier  (DSB-SC) 

Amplitude  Modulation  (AM) 

b) VHF  band  (30  -  300  MHz):  (with  25  KHz  channel  spacing) 

Frequency  Modulation  (FM)  :  deviation  5.0  to  5.6  KHz 
Non-Coherent  Frequency-Shift  Keying  (XC-FSK)  :  16  Kbps 
and  19.2  Kbps  for  Frequency  Hopping 

there  is  a  potential  of  [(27M/3K)+(270M/25K)]  -  19800  distinct  modulated  signals.  It 
is  therefore  permited  to  assume  a  maximum  bandwidth  per  signal  of  at  most  30  KHz. 
Hence,  in  order  to  isolate  a  SOI  from  the  other  signals,  a  means  of  digitally  filtering  and 
mixing  the  signal  to  baseband  is  required.  In  addition,  the  signal  must  be  decimated  to 
provide  a  signal  suitable  for  further  processing  by  conventional  DSP  devices  such  as  a 
TMS320C30.  This  last  operation  must  give  a  new  sampling  frequency  b1(H,,v,  being 
approximately  twice  the  bandwidth  of  any  expected  SOI  bandwidth  tsav  30  KHz). 
Therefore,  fs(now)  will  be  something  around  60  KHz,  an  attainable  speed  tor  conventional 
DSP  chips. 
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2.0  DIGITAL  FILTER  CONSIDERATIONS 


2.1  CHOICE  OF  FILTER 

Because  it  is  desired  to  keep  the  phase  information  intact,  a  linear  phase  filter  has 
to  be  used.  Therefore,  the  Finite  Impulse  Response  -  FIR  -  (non-recursive)  filter  is  more 
appropriate  than  an  Infinite  Impulse  Response  -  IIR  -  (recursive)  filter.  Also,  it  offers  the 
advantage  of  being  always  stable  whereas  it  is  not  the  case  for  an  IIR  filter. 
Furthermore,  the  quantization  noise  is  not  additive  (it  is  always  the  same  for  every 
output).  Some  FIR  filter  characteristics  will  be  discussed  below. 


First,  the  sharpness  of  the  filter  is  proportional  to  the  filter  length.  The  more  taps 
the  filter  has,  the  more  sharper  band  edge  characteristic  it  will  have.  Second,  the  filter 
window  type  dictates  the  kind  of  behaviour  the  filter  will  have  inside  the  passband 
region  as  well  as  in  the  stopband  region.  The  Kaiser  window  (KW)  is  one  of  the  best 
types  of  window  because  it  provides  a  parameter  called  beta  ((5)  which  allows  a 
compromise  between  transition  band  requirement  and  stop  band  attenuation.  The  KW 
impulse  response  is  expressed  as  follows: 


w(n) 


P 


'M-\V  f 


n  - 


M- 1 


V 


M- 1 


P 


,  0  <  «  <  M  ■ 


(1) 


where  M  is  the  filter  length  (or  number  of  taps),  n  the  sampling  instant  (an  integer)  and 
I„  is  the  modified  Bessel  function  of  order  zero.  Equations  that  relate  (3  to  the  stopband 
attenuation  A  are 

(3-0  ,  A<2\dB 

-  0.5842  [A -2 1)° 4  +  0.07886 (a -2 1)  ,  2\ dB  <  A  <  50 ciB  <2) 

-0.1 102 (4  -8.7)  ,A>50dB. 

From  Eq.2,  it  is  easy  to  compute  (3  that  corresponds  to  a  stopband  attenuation  of  let's  sav 
80  dB.  One  will  find  (3  =  7.64.  Fig.l  pictures  the  shape  of  the  impulse  response  of  the 
KW  given  in  Eq.l  for  different  (3  (4,  6,  7.64)  and  M  =  255. 

It  was  found  that  a  filter  of  length  M  =  2o5  using  a  KW  with  (3  =  5.66  provides 
adequate  characteristics  to  the  filter:  60  dB  stopband  attenuation  and  a  nice  sharp  band 
edge  as  seen  in  Fig. 2.  Note  that  if  P  is  reduced,  the  transition  width  goes  down  but  the 
stopband  attenuation  is  also  diminished. 

The  KW  gives  the  filter  designer  the  flexibility  that  other  tvpes  of  windows  like 
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Rectangular,  Bartlett,  Hanning,  Hamming  and  Blackman  do  not  have.  All  ot'  those  have 
a  fixed  peak  side  lobe  and  a  fixed  equation  available  for  the  transition  width  of  main 
lobe  as  shown  in  Table  I. 


Recall  that  the  transition  width  and  the  peak  side  lobe  tor  the  KYV  filter  design 
do  not  only  depend  on  the  filter  length  M  but  also  depends  on  the  parameter  (3. 
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Table  I:  Frequency-Domain  Characteristics  of  Some  Window  Functions 


Type  of  Window 

Transition  Width 
of  Main  Lobe 

Peak  Side  Lobe 
(dB) 

Rectangular 

4k/M 

-13 

Bartlett 

8k/M 

-27 

Hanning 

8k/M 

-32 

Hamming 

8k/M 

-43 

Blackman 

12k/M 

-58 

2.2  FILTER  LIMITATIONS 

Before  going  further,  a  more  convenient  normalized  frequency  has  to  be  defined. 
In  the  digital  domain,  it  is  always  helpful  to  use  a  normalized  digital  frequency  denoted 
9  which  is  equal  to  (2k- f/Q.  From  now,  it  will  be  easy  to  go  from  f  to  0  or  vice-versa 
whenever  needed  (depending  on  which  form  is  best  suited  for  a  particular  problem). 

The  major  limitations  of  the  filter  are  caused  by  the  transition  width  and  the 
stopband  attenuation.  The  transition  width  limits  how  narrow  the  width  of  the 
passband  can  be.  Also,  with  a  KW  filter,  adjusting  p  such  that  the  stopband  attenuation 
required  is  met  affects  the  transition  width.  It  is  always  possible  to  diminish  the 
transition  width  by  increasing  the  number  of  taps  M,  but  in  real-time  DSP,  it  is,  for  a 
number  of  reasons,  really  difficult  to  do  that  kind  of  operation  after  a  certain  value  of 
M  has  been  reached.  Therefore,  M  has  to  be  as  high  as  possible  but  it  has  a  certain 
practical  hardware  limitation. 

A  filter  bandwidth  limit  can  be  found  assuming  a  KW  filter  with  255  taps  and  60 
dB  stopband  attenuation.  For  the  discussion,  the  transition  width  is  defined  as  follows: 

O<0  <0,  — >4  < 0. 1  dB 

02<0<k  > 60 dB  (3) 


assuming  that  a  low-pass  filter  is  used.  The  width  of  transition  is  easily  found 
experimentally;  for  example,  a  KW  low-pass  filter  having  the  above  characteristics  has 
been  designed  with  a  3  dB  cut-off  frequency  equal  to  0.5ji  (=f„/4).  It  gave  the  following 
results: 
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0,  -0.490k  ,  0  -0.517  k  ,  A0 

1  ’2  ’  transition 


-  0.027  k  . 


(4) 


This  last  result  is  of  great  importance  and  was  found  to  be  independant  of  the  filter  cut¬ 
off  frequency.  It  means  that  the  filter  must  add  0.054k  (which  is  2  times  the  transition 
width)  to  the  bandwidth  required  for  the  SOI  (something  around  30  KHz).  This  also 
means  that  the  minimum  filter  bandwidth  m  0.054k. 

2.3  FILTER  SPECIFICATIONS 

The  specifications  for  all  filters  used  in  next  sections  are  described  by  the 
following  lines: 

—  Kaiser  Window  Filter 
P  =  5.66 

Transition  width  =  0.027k 
Stopband  attenuation  =  60  dB 
Filter  length:  M=255 


It  represents  a  good  compromise  between  stopband  attenuation,  transition 
bandwidth  and  the  filter  length  M.  Note  that  the  transition  width  as  defined  in  Eq  3  is 
not  a  usual  way  to  define  it  since  for  filter  design  purpose,  the  cut-off  fr-y^ncy  is 
defined  as  being  the  3  dB  attenuation  point.  Therefore,  it  was  found  experimentally  that 
for  this  particular  KW  filter,  going  from  0.1  dB  to  3  dB  attenuation  point  corresponds  to 
a  width  of  0.0098k. 

From  now  on,  no  matter  what  filter  type  will  be  used  (low-pass,  highpass, 
bandpass,  etc...),  the  above  specifications  will  always  hold. 


3.0  OVERALL  DESIGN  STRATEGY 

3.1  ENERGY  DETECTION 

A  wideband  digital  receiver  located  at  the  front  end  of  an  intercept  system  is 
assumed  to  be  capable  of  looking  at  frequencies  from  3  MHz  up  to  at  least  300  MHz 
with  4  MHz  of  bandwidth  at  a  time.  Hence,  the  first  operation  after  the  reception  is 
called  Energy  Detection  (ED).  It  has  to  determine  if  and  where  a  SO!  is  present  inside 
the  4  MHz  bandwidth.  If  yes,  then  it  detects  an  approximative  center  frequ  mcv.  At  the 
present  time,  the  ED  is  seen  as  a  system  having  a  frequency  resolution  of  1  KHz.  The 
Mixing,  Filtering  and  Decimation  System  (MFDS)  receives  'he  approximate  center 
frequency  from  the  ED  system  and  it  then  begins  to  isolate  the  SOI. 
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3.2  MIXING,  FILTERING  AND  DECIMATION  SYSTEM 


Following  the  discussion  from  Section  2,  it  is  known  that  filtering  operations 
restrict  the  design  of  the  MFDS.  For  a  sampling  frequency  of  fs  =  8  MHz,  the  minimum 
filter  bandwidth  is  [(0.0547t-8M)/2jt)]  =  216  KHz.  Adding  the  SOI  bandwidth  (^ssumed 
to  be  30  KHz)  to  this  transition  bandwidth  makes  a  total  filtering  bandwidth  of  246  KHz. 
It  is  therefore  obvious  that  there  is  no  way  to  go  directly  from  a  sampling  frequency  of 
8  MHz  to  ~60  KHz  in  only  one  decimation.  It  will  have  to  be  done  in  two  major  similar 
operations. 

Briefly,  the  steps  taken  through  this  process  will  be  the  following  ones: 

1-  ED  gives  fcl  (center  frequency)  of  a  SOI. 

2-  Complex-mix  from  fcl  to  f^  =  500/4  KHz,  followed  by  filtering  (bandwidth 
of  complex  bandpass  filter  is  BW1  =  246  KHz  ’). 

3-  Decimate  by  16  to  go  from  fsl  =  8  MHz  to  fs2  =  500  KHz,  the  sampling 
frequencies. 

4-  Complex-mix  from  fc2  to  f^  =  83.33/4  KHz  and  filter  (B\V2  =  43.5  KHz  2). 

5-  Decimate  by  6  to  go  from  fs2  =  500  KHz  to  f^  =  83.33  KHz. 

The  last  sampling  frequency,  f^  =  83.33  KHz,  is  the  minimum  sampling  frequency 
that  is  feasible  to  reach  in  order  to  insure  that  a  stopband  attenuation  of  60  dB  is 
respected  through  all  decimation  processes  and  the  30  KHz  bandwidth  for  the  SOI  has 
less  than  0.1  dB  attenuation.  Those  steps  guarantee  a  very  good  protection  against 
possible  aliasing  caused  by  the  mixing  and  decimation  operations.  Note  that  the  center 
frequency  is  redefined  at  steps  2-  and  4-  to  be  equal  to  the  new  sampling  frequency 
divided  by  4  (ie.,  fs/4).  The  next  section  will  focus  on  the  details  of  each  individual  part 
of  this  process. 


4.0  DETAILED  SYSTEM  DESCRIPTION 

4.1  COMPLEX  FILTER  DESIGN 

The  filter  has  to  be  designed  first  as  a  low-pass  filter  with  a  rectangular  window 
transfer  function  having  the  impulse  response  given  by 


This  bandwidth  is  the  total  tiller  bandwidth  where  the  attenuation  is  less  than  the  stopband  attenuation  (oOdR) 
Note  that  the  (liter  bandwidth  used  for  the  design  is  |30K-H0.009X;r-2-XM,'2"i|  =  108.4  KM/  This  result 
corresponds  to  the  filter  bandwidth  defined  at  the  3  dB  attenuation  points. 

BW2  =  [30K+(0.054jf500K)/2jt|  =  43.5  KH/.  As  in  note  1.  the  3  dB  tiller  bandwidth  is 
[30K+«).00987i-2-500K/2;t)|  =  34.9  KH/. 
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,  j/  M  is  selected  to  be  odd. 
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It  is  already  known  that  the  filter  length,  M,  will  be  odd  to  insure  a  symmetric 
filter  spectrum  as  opposed  to  an  anti-symmetric  one  when  the  even  length  is  used. 
Therefore  Eq.6  will  always  hold.  The  relationships  between  fc,  ooc  and  0„vv  are  shown 
by  the  following  equations: 


to 
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2nfc 

~r 
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BW 


20) 


(7) 


where  0BW  is  the  total  filter  bandwidth. 

Those  three  last  equations  define  a  Rectangular  window  filter  with  cut-off 
frequency  coc  being  in  radians  and  chosen  between  0  and  n.  Mow,  the  filter  impulse 
response  h(n)  has  to  be  multiplied  by  the  Kaiser  Window  impulse  response  w(n)  given 
in  Eq.l.  (J  has  to  be  set  at  5.66  to  provide  specifications  as  discussed  earlier.  The  filter 
impulse  response  is  then  expressed  as 

hu>(n)  *  h{n)w{n).  (8) 


In  the  filtering  operation,  the  filter  required  is  a  complex  bandpass  filter.  The  way  of 
working  with  such  a  filter  is  to  separate  the  real  and  imaginary  part  of  a  complex  low- 
pass  filter.  The  bandpass  complex  filter  is  created  as  follows: 


hc{n)  -  hLP(n)eJ*’n  -  hr(n)  +jhi(n) 

hr{n)  -  hLP(n)cos(dcn)  -> Hr(Q )  -  -  0f)  +  HJ?  ♦  0c) 

ht(n)  -  hLP{ri)  sin(0c  n)  -4^(0)  -  J_[//^0  -  0c)  -  Htn(Q  +  0  ) 


(9) 


where  -4  stands  for  the  Fourier  Transform  operation. 


The  parameter  0C  is  the  bandpass  filter  center  frequency.  Note  that  the  power 
spectrum  density  (PSD)  of  Hr(0)  and  H,(0)  are  identical  but  their  phases  are  not.  See 
Figs.  3,  4  and  5  for  visual  representations  of  respectively  Hr(0)  and  A.Hr(0),  H,( 0)  and 
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4.H,(0),  Hc(0)  and  AHc(0). 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8 

Frequency  (xPl  radians) 


Figure  3:  Power  Spectrum  And  Phase  Plot  Of  hr(n) 

In  this  example,  0C  has  been  chosen  to  be  0.5tt  and  the  filter  bandwidth  BW  = 
0.4tc.  On  the  graphs,  the  value  1  on  the  x  axis  coincides  with  fs/2.  The  part  going  from 
0  to  1  is  the  positive  frequency  side  while  the  part  from  2  to  1  is  the  negative  frequency 
side  (2  corresponding  to  the  0  frequency  and  1,  the  -fs/2  frequency).  Recall  that  in  the 
digital  domain,  the  spectrum  is  always  periodic  with  a  period  being  equal  to  fs  (or  2k 
in  0  domain). 
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frequency  (xPl  radians) 


0  0.2  0.4  0.6  0.8  1  1.2 

frequency  (xPl  radians) 


Figure  4:  Power  Spectrum  And  Phase  Plot  Of  h,(n) 


Frequency  (xPl  radians) 

Figure  5:  Power  Spectrum  And  Phase  Plot  Of  hc(n) 
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4.2  COMPLEX  MIXING 


The  complex  spectral  shifting  block  diagram  is  shown  in  Fig.6. 


The  input,  x(n),  is  the  signal  to  be  mixed-down.  The  output  signals  are  in-phase 
and  quadrature  components  of  the  frequency  shifted  signal,  ie.,  respectively  i(n)  and 
q(n).  As  it  can  be  seen  from  Fig.6,  the  input  signal  x(n)  is  mixed-down  to  0C  -  05hlfl  = 
9c(new),  a  new  center  frequency,  assuming  that  x(n)  has  a  certain  center  frequency,  0C.  The 
filters  hr(n)  and  h;(n)  are  the  real  and  imaginary  parts  of  the  complex  filter  as  defined 
in  Eq.9  but  with  0C  now  being  changed  to  0c(new).  Make  note  that  the  same  filter  is  used 
twice  within  this  operation.  The  bandwidth  of  the  input  signal  is  0„w.  The  filter  cut-off 
frequency  that  has  to  be  used  for  Eqs.  5  and  6  is  based  on  Eq.7  and  is 


0 


0)  - 


BW 

“ 


(10) 
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The  general  theory  under  this  complex  mixing  scheme  can  be  found  in  Ref.[l], 
Section  2.2.  The  next  set  of  equations  (Eq.ll)  briefly  presents  all  the  steps  one  has  to  go 
through  to  arrive  at  the  complex  mixing  block  diagram  shown  in  Fig. 6.  Notice  that  the 
symbol  *  denotes  a  convolution  operation. 

hc(n)  -  hLP(n)eJ*«-'n  -  hr(n)  +jh(n) 
i\n)  -  2x(n)cos(Qski/in) 
q’{n)  -  -24«)sin(e^n) 

x\n)  -  i\n)  +  jq'(n)  (11) 

y(n)  -  x'(n)  *  hc(n)  -  [/'(«)  +jq\n)]  *  \hr(n)  +J/7i,(/i)] 


i(n)  -  9?{y(/i)}  - 

i'{n)  *  hr(n) 

- 

q\n)  *  ht(n) 

q(n)  -  3{y(«)}  - 

i'(n)  *  ht(n) 

+ 

q\n)*hr(n) 

As  one  may  observe,  if  q(n)  is  not  wanted,  the  number  of  parts  could  be  reduced 
to  half.  That  will  make  considerable  hardware  cost  savings  when  implementing  the 
complete  MFDS  but  that  will  be  discussed  later. 

4.3  DECIMATION 

A  decimation  is  a  fairly  simple  operation.  What  seems  to  be  more  complex 
constitutes  any  process  preceding  it.  The  conditions  have  to  be  such  that  no  aliasing  can 
arise  under  the  decimation.  Therefore,  a  special  care  should  be  taken  during  the  mixing 
and  filtering  processes  as  explained  in  Section  3.  Also,  the  decimation  factors  have  to 
be  scrupulously  chosen.  In  Section  3.2,  those  factors  were  determined  to  be  16  and  6. 
Those  values  are  selected  to  maximize  the  sampling  frequency  reduction  without 
affecting  the  SOI  bandwidth. 

To  decimate  the  signal  by  16  or  6,  the  samples  have  to  be  picked  1  every  16  or  1 
every  6  samples  of  data  respectively. 


5.0  COMPLETE  SYSTEM  IMPLEMENTATION 

5.1  BLOCK  DIAGRAM 

The  block  diagram  of  the  MFDS  is  shown  in  Fig. 7.  It  basically  reproduces  what 
has  been  defined  in  Section  3.2. 

To  compute  the  shifting  frequencies,  0shllll  and  0shll,:,  Eq.12  is  used,  ie.. 
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Figure  7:  Block  Diagram  Of  Complete  MDFS 
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5.2  DISCUSSION 


The  upper  part  of  the  block  diagram  shown  in  Fig. 7  represents  the  first 
decimation  process  changing  a  sampling  frequency  of  8  MHz  to  500  KHz.  This  part  onlv 
needs  to  produce  the  in-phase  component  since  the  quadrature  component  is  not 
required  as  the  input  to  the  second  decimation  process  (lower  part  of  Fig. 7).  As  a  result, 
half  the  filters  may  be  dropped  from  the  first  decimation.  The  second  decimation 
changes  the  500  KHz  sampling  frequency  to  83.33  KHz  and  provide  outputs  that  are  the 
in-phase  and  quadrature  components,  i(n)  and  q(n).  Let's  look  more  in  detail  at  the 
hardware  speed  requirements  for  each  part. 
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The  two  first  mixers  need  to  compute  one  multiplication  in  [1/8M]  =  125  ns. 
Most  common  DSP  processors  can  do  that  job  easily.  Next,  the  two  filters  have  to 
calculate  the  following: 

M- 1 

FilterOut  -  Y'  bFilterln  (13) 

m-0 

where  bm  are  the  filter  weights.  As  seen  from  Fig. 7,  it  looks  like  the  filter  outputs  data 
at  the  same  rate  as  it  comes  in  ie.,  8  MHz.  But,  the  filtering  and  decimating  processes 
can  be  incorporated  such  that  the  filter  output  rate  could  be  fs2,  500  KHz.  To  do  so,  the 
filter  computes  an  output  data  at  {  n=0,16,32,...  }.  Therefore,  it  must  have  a  queue  to 
keep  buffered  16  data  before  processing  them.  Hence,  the  equation  will  be: 

254 

FilterOut  -'V'  bFilterln  ,  h-0,  16,32,...  d1*) 

n  171  n—rrl 

m-0 

Eq.14  involves  M  (=255)  multiplications  and  M  additions.  Therefore,  2M 
operations  must  be  performed  so  that  it  takes  [1  / (500K-2-255)]  =  3.92  ns  /  operation 
for  each  filter.  At  best,  assuming  a  processor  able  to  do  add /multiply  in  the  same 
instruction,  then  it  would  take  [3.92ns-2]  =  7.84  ns  /  instruction.  Actual  common  DSP 
chips  are  not  able  to  process  signals  at  that  speed  (eg.,  TMS320C30  has  60  ns  single 
instruction  cycle).  Some  state-of-the-art  processor  boards  like  the  Intel  iVVarp  processor 
boards  (not  considered  here  as  common  DSP  boards)  could  be  an  alternative  but  they 
are  very  expensive.  Nevertheless,  they  are  really  fast  and  very  flexible.  For  real-time 
intensive  computation  as  for  DSP  applications,  those  boards  are  really  well  designed. 
Each  iWarp  processor  is  composed  of  a  computation  agent  and  a  communication  agent. 
These  agents  are  working  simultaneously  such  that  transfer  delays  between  two 
processors  are  minimized.  Thus,  making  a  very  powerful  parallel  computer  with 
iWarp's  is  something  possible.  Using  TM5  family  of  processors  will  never  reach  the 
degree  of  parallelism  one  can  attain  using  iWarp's. 

The  500  KHz  hjfn)  output  data  is  substracted  from  the  output  of  the  hr(n)  filter. 
It  requires  a  maximum  instruction  cycle  of  [1/500K]  =  2  ps.  Such  a  speed  is  easy  to 
reach  with  a  TMS320C30.  The  same  comment  applies  for  the  two  mixers  in  the  second 
decimation  system.  The  filtering  can  be  done  along  with  the  decimation  as  in  the  first 
decimation  system.  It  is  possible  to  compute  the  required  processing  speed  ie., 
[1  /83.33K*255]  =  47  ns  /  instruction.  Again,  use  of  iWarp  processors  may  be  best  suited 
in  this  case.  The  two  adders /subtractors  at  the  end  must  work  at  [1/83. 33K]  =  12  ps  / 
addition. 

At  the  output  of  the  system,  i(n)  and  q(n)  are  provided  at  a  rate  of  83.33  KHz 
allowing  a  signal  bandwidth  equal  to  41.67  KHz.  a\  visualization  of  the  whole  MFDS 
process  is  shown  in  Fig.8. 
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Each  numbered  arrow  symbolizes  a  particular  part  of  the  MFDS  operation. 
Arrow  #1  represents  the  first  down  mixing  and  filtering,  #2,  the  first  decimation,  #3,  the 
second  down  mixing  and  filtering  and  finally,  #4  denotes  the  last  decimation.  The  last 
graph  in  the  bottom  right  corner  pictures  the  PSD  of  i(n)  or  q(n).  The  limit  of  this 
system  resides  in  a  signal  band  limited  to  41.67  KHz.  It  may  therefore  be  possible  to 
encounter  more  than  one  SOI  inside  this  bandwidth  according  to  the  possible  types  of 
modulation  previously  described.  Additional  filterings  will  have  to  be  made  to  isolate 
only  one  SOI.  Nevertheless,  the  MFDS  is  fundamentally  designed  to  decrease  the  data 
sampling  frequency  such  that  a  DSP  chip  like  a  TMS320C30  will  be  able  to  continue  anv 
required  further  processes.  Anyhow,  the  theory  behind  the  MFDS  may  serve  as  a 
reference  for  complex  down  mixing  and  decimation  operations.  Note  from  Fig.8  that  the 
signal  on  the  fifth  graph  is  seen  as  a  signal  having  an  unknown  bandwidth  frequency 
BW3  and  also,  fc3  constitutes  an  approximate  center  frequency.  All  post-processing  is  not 
considered  here  because  it  is  out  of  scope  of  this  technical  memorandum. 
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6.0  SIMULATIONS  AND  RESULTS 


6.1  SIMULATION  SETUP 

The  present  simulation  is  intended  to  demonstrate  the  validity  of  the  theory  used 
under  MFDS.  For  this  matter,  only  one  complex  mixing,  filtering  and  decimation  will 
be  performed  instead  of  two  as  in  MFDS  since  the  operations  are  the  same.  Also,  the 
values  will  not  be  identical  to  what  was  previously  discussed.  The  following  describes 
the  simulation  conditions: 

The  input  signal,  xn,  is  a  block  of  data  of  length  16384  having  a  sampling 
frequency  of  fsl  =  240  KHz.  It  is  composed  of  three  distinct  signals:  one  is 
a  BPSK  signal  at  a  center  frequency  fc,  =  90  KHz  with  a  bandwidth  of  6 
KHz,  the  second  one  is  an  ASK  signal  at  fc  =  50  KHz  with  a  bandwidth  of 
15  KHz  and  finally  the  last  one  is  another  BPSK  signal  at  fc  =  110  KHz 
with  a  bandwidth  of  10  KHz.  The  modulation  envelopes  were  generated 
from  digitized  samples  of  an  actual  voice. 

The  first  BPSK  signal  located  at  fcl  =  90  KHz  and  BW1  =  6  KHz  will  be 
considered  as  the  SOI  to  isolate. 

The  objective  is  to  take  that  BPSK  signal  at  a  240  KHz  sampling  rate,  mix 
it  down  to  a  new  center  frequency  fc2  =  7.5  KHz  and  decimate  it  to  get  a 
new  sampling  frequency  fs2  =  [240K/8]  =  30  KHz  (decimation  bv  8). 

With  the  same  data  used  to  generate  the  SOI,  a  BPSK  signal  with  a 
sampling  frequency  of  30  KHz,  a  bandwidth  of  6  KHz  and  a  center 
frequency  of  7.5  KHz  will  be  built.  Through  the  simulation,  this  signal  will 
be  called  B3  and  will  serve  as  a  reference  to  compare  to  the  system  output, 
Yd.  Y  will  stand  for  the  down  converted  signal  before  decimation.  To 
compare  Y  to  its  expected  result,  an  other  BPSK  signal  called  B4  will  be 
generated  at  fc2  =  7.5  KHz  with  a  6  KHz  bandwidth  and  a  sampling 
frequency  of  240  KHz. 

For  filter  design  purposes,  the  3  dB  cut-off  frequency  is  equal  to 
[6K+(0.0098'7t-240K/jt)]  =  8.352  KHz.  Also,  the  total  bandwidth  for  60  dB 
attenuation  will  be  [6K+(0.0547tT20K/rc)]  =  12.48  KHz  which  has  to  be 
lower  than  fs/2  (15  KHz). 

6.2  SIMULATION  ALGORITHM 

The  present  section  describes  a  complete  algorithm  illustrating  the  program  steps. 
It  is  proposed  as  a  guideline  for  any  software/hardware  implementation  of  a  complex 
mixing,  filtering  and  decimation  system. 
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Step  1: 


-  Inputs: 


x  =>  block  of  data  of  length  N1 =16384 

sf  =>  sampling  frequency 

cf  =>  center  frequency  of  the  SOI 

bwf  =>  3  dB  bandwidth  allowed  for  filter  according  to  SOI 
bandwidth 

ncf  =>  new  center  frequency 
dr  =>  decimation  ratio 
M  =>  filter  length 

(3  =>  beta  parameter  for  KW  =  5.66  for  60  dB  attenuation 

Step  2:  -  Complex  Mixing: 

0shift  =  2rr(cf-nc0/sf 
0ncf  =  2;r(ncf)/sf 
FOR  n=0  to  N 1-1,  DO 

xr(n)  =  2x[n}cos(Qshifn) 

x,(n)  -  -2x[n]  sin (Qsltifn) 

END 

Step  3:  -  Filter  Design: 

toc  =  7r(bwf)/sf 

FOR  n=0  to  [((M-l)/2)-l],  DO 


h[n) 


w[«]  - 

hLP{n) 
hr{n)  - 

h'(n)  -  hLp(n)sm(Qnr/n) 

END 
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4 


2 

0) 

hLP(P)  -  — 
n 

hf(P)  ’  hLP^COS(Qncfn) 

FOR  n=[((M-l)/2)+l]  to  M-l,  DO 

/ir(«)  -  hr(Nl~n) 
ht(n)  -  h^Nl -n) 

END 


Step  4:  -  Filtering: 


j  -  cconv[hr ,  jcr)  - 
q  -  cconv(hi ,  jtr)  + 


cconv(ht ,  .vj 
cconv(hr ,  x^j 


Step  5:  -  Decimation: 


FOR  n=0  to  [(Nl/dr)-l],  DO 

idec(n)  -  i(drxn) 
qdec(n)  -  q(drxn) 

END 

Step  6:  -  Outputs: 

i  =>  in-phase  output  component  before  decimation 
q  =>  quadrature  output  component  before  decimation 
idee  =>  in-phase  output  component  after  decimation 
qdec  =>  quadrature  output  component  after  decimation 

Note  that  in  Step  4,  the  function  "cconv”  performs  a  circular  convolution  between 
the  two  arguments  in  parenthesis.  The  output  of  a  circular  convolution  has  the  same 
length  of  the  argument  that  has  the  greater  length.  The  formula  for  a  circular 
convolution  is: 


M- 1 

y(n)  -  ^  h(m)x[MOD(n-m,N)\. 

m-0 
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for  n=0,l,...,N-l  where  N  is  the  length  of  x  and  M,  the  length  of  h.  The  MOD  function 
finds  the  modulo  N  value  of  the  first  argument.  The  circular  convolution  is  used  here 
because  the  operation  involves  a  block  of  data  manipulations  instead  of  a  continous  real¬ 
time  data  stream. 

6.3  RESULTS  OF  SIMULATIONS 

In  this  section,  the  power  spectral  density  (PSD)  of  the  signals  will  be  compared 
(in  log  scale).  On  the  graphs,  the  following  symbols  will  be  used: 

Y  =>  PSD  of  i(n),  the  in-phase  component  before  decimation. 

X  =>  PSD  of  x(n),  the  input  signal. 

Yd  =>  PSD  of  idec(n),  the  in-phase  component  after  decimation. 

B3  =>  PSD  of  directly  generated  SOI  signal  at  a  30  KHz  rate. 

B4  =>  PSD  of  directly  generated  SOI  signal  at  a  240  KHz  rate. 


As  seen  from  Fig. 9,  X  shows  clearly  the  three  distinct  signals  between  0  and  120 
KHz.  Also,  because  the  signal  is  real,  the  spectral  svmmetrv  characteristic  can  easily  be 
observed.  The  SOI  is  in  the  middle  of  the  three  first  left  signals  (ie  ,  at  f,  =  ^0  K1  lz)  and 
its  symmetrical  image  is  seen  at  f.  =  150  KHz  or  [150K-240K]  =  -*->0  KHz. 

The  result  of  the  complex  mixing  and  filtering  operation,  depicted  in  Fig.  10  as  V, 
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shows  only  one  non-attenuated  signal  which  is  the  mixed  and  filtered  SOI  at  7.5  KHz. 
This  signal  can  also  be  seen  at  -7.5  KHz  or  [240K-7.5K]  =  232.5  KHz.  All  the  other 
signals  are  attenuated  by  at  leat  60  dB.  A  comparison  between  Y  and  B4  is  made  in 
Fig.ll  but  it  is  hard  to  make  any  judgement  on  the  result  since  only  6  KHz  is  a 
bandwidth  of  interest  (BOI).  Therefore,  a  zoom  of  this  particular  band  (4500  -  10500  Hz) 
has  been  illustrated  in  Fig.12.  In  those  two  last  figures,  Y  is  seen  as  a  solid  line  while 
B4  as  a  dotted  line.  From  Fig.12,  only  very  tiny  discrepancies  are  noticed.  The  absolute 
value  of  the  error  between  Y  and  B4  within  the  BOI  is  shown  in  Fig.  13.  The  mean  of  this 
error  is  found  to  be  ~0.4  dB. 
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After  the  decimation  has  been  performed.  Yd  and  B3  are  compared  in  Fig. 14.  A  % 

closer  look  at  it  is  taken  in  Fig.15  (0  -  fs/2).  Finally,  the  same  frequency  band  as  in 
Fig.12  is  used  to  compare  Yd  and  B3  in  Fig. 16.  An  increased  error  is  observed  from 
there.  Now,  the  mean  of  the  error  between  Yd  and  B3  is  equal  to  1.1487  dB.  Again,  an  * 

error  graph  within  the  BOI  is  drawn  in  Fig.  17. 
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Figure  17:  Error  Function  Between  Yd  And  B3  Within  The  BOI 


An  important  observation  has  to  be  made  from  the  above  results.  The  value  of 
the  error  is  less  meaningful  than  it  appears.  As  it  should  be  noted  from  Figs.  12  and  15, 
the  errors  appear  where  the  valleys  are  and  not  at  the  peaks.  Therefore,  errors 
encountered  through  the  simulation  should  be  considered  negligible.  Given  this  fact,  the 
validity  of  the  theory  used  under  MFDS  has  been  confirmed.  The  observed  error  is 
mainly  due  to  the  aliasing  coming  from  the  60  dB  attenuated  signals. 


7.0  CONCLUSION 

This  work  had  two  objectives.  First,  to  reduce  the  sampling  frequency  given  bv 
a  digital  receiver  in  the  ACES  system  front  end  to  an  attainable  speed  for  a  T.V1S320C30 
chip.  Second,  isolate  a  SOI  without  altering  it.  These  two  objectives  were  presented 
conceptually  in  Section  6  through  computer  simulations.  Prior  to  this  section,  a  lot  of 
work  was  done  to  determine  the  filter  specifications,  define  a  complex  mixing  scheme 
and  incorporate  a  decimation  operation  to  complete  an  overall  system  design. 

A  set  of  parameters  was  assumed  for  the  ACES  system  in  order  to  make  the 
problem  as  practical  as  possible  at  this  stage  (eg.,  a  receiver  with  8  MHz  sampling  rate). 
One  can  also  refer  to  Fig. 8.  Speed  requirements  were  discussed  in  Section  3  and  from 
there,  it  is  anticipated  that  some  kind  of  really  fast  parallel  processing  like  iW'arp  boards 
will  be  needed.  Customized  hardware  could  also  be  used  as  an  alternative  but  was  not 
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mentioned  earlier  because  that  approach  would  be  too  restrictive. 

Briefly,  the  MFDS  can  be  seen  as  two  non-coherent  mixings.  Unfortunately,  the 
SOI  is  not  completely  filtered  at  the  end  and  it  is  therefore  envisioned  that  further 
processing  will  be  required  before  studying  this  particular  signal.  Consequently,  a 
significant  amount  of  work  will  be  carried  on  to  develop  a  system  called  a  Digital 
Coherent  Adaptive  Down  Converter.  The  output  of  this  system  will  have  to  be  the  SOI 
only,  completely  filtered,  with  a  center  frequency  of  0  and  possibly  a  reduced  sampling 
rate. 
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